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The influence functional method of Feynman and Vernon is used to obtain a quantum master 
equation for a Brownian system subjected to a Levy stable random force. The corresponding classical 
transport equations for the Wigner function are then derived, both in the limit of weak and strong 
friction. These are fractional extensions of the Klein-Kramers and the Smoluchowski equations. It 
is shown that the fractional character acquired by the position in the Smoluchowski equation follows 
from the fractional character of the momentum in the Klein-Kramers equation. Connections among 
fractional transport equations recently proposed are clarified. 
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In the theory of Brownian motion one is interested in the time evolution of a system coupled to a large environment. 
The effect of the coupling is modeled by a stochastic force F{t) with a given probability density P[F(t)]. The dynamics 
of a Brownian particle of mass M in the presence of an external potential U(x) is then described by the Langevin 
equation 

Mx(t) + jx(t) + U'(x) = f (i) , (1) 

where F(t) has been divided into a mean force proportional to the velocity the friction force — 7x(t), plus a fluctuating 
part In the usual treatment of Brownian motion it is assumed that the random force is Gaussian distributed 
with variance = 2D 5[t — t') where D = jkT is the diffusion coefficient, and the Langevin equation is shown 

to be fully equivalent to a phase-space equation — the Klein-Kramers equation. In the limit of strong friction, the 
inertial term in the Langevin equation can be neglected and the Klein-Kramers equation reduces to the Smoluchowski 
equation. However, it has become clear in recent years that many processes in nature, such as anomalous diffusion 
(for a review see |^-^]), cannot be described by ordinary (Gaussian) Brownian motion. A case in point is the so-called 
Levy flight with a stochastic force distributed according to Levy stable statistics and which has been introduced in 
connection with super-diffusion [^|^] . In this Letter we consider the generalization of transport equations to describe 
Levy Brownian motion. This question has already been addressed in the past by using various methods in 
particular the CTRW formalism |^jn|. However, all these approaches were limited to coordinate space only. Here we 
present the first derivation of the Klein-Kramers equation for a Levy stable process. We consider both the case of a 
symmetric and asymmetric probability distribution. As our main tool, we employ the influence functional formalism 
developed by Feynman and Vernon [ |ll|Jl2| . 

If initially system and environment are not correlated then, according to Feynman and Vernon 0Jl|, the density 
operator of the system at time t can be written in coordinate representation as 

p(x,x',i) — / J-[x,x') exp — {s[x}-S[x']jp(x o ,x o ',0) Vx(t)Vx'(t)dx dx ' . (2) 



Here the entire information on the coupling to the environment is contained in the influence functional F[x, x'] 
where is the characteristic functional of the probability density P[F(t)], 



$[&(*)] = J exp J F(t)k(t)dt} P[F{t)]VF{t) . (3) 

If F(t) is Gaussian distributed with mean (F(t)) = f(t) and variance {[F(t) - j(t)][F(t') - 7(f)]) = D(t - t'), the 
characteristic functional is given by 

$ G [fc(*)] = exp ji J j(t)k(t)dt -~f j D(t - t')k(t)k(t')dtdt'} . (4) 
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If we further assume that the friction force is proportional to the velocity of the system, 7(f) = —j[x(t) + x'(t)]/2, 
and that the variance is delta correlated in time, D(t — t') = 2D S(t — t'), then the influence functional can be written 
in the form 



6XP ^ 



1 [x(t) - x'(t)][x(t) + x'(t)} + i j [x(t) - x'(t)] 2 }dt . 



(•5) 



By means of a small time expansion, Eq. (g) can be transformed into a differential equation for the density operator. 
Using the influence functional (|^), this results in the master equation 



ih 



dp(x, x' , t) 
di 



H[x) - if (x') + -^{x- x')(p x -p x >)- 1 j (x - x') 2 ] p(x, x', t) , 



(6) 



where H(x) — p x 2 /2M + U(x) is the Hamiltonian of the system. We recognize in Eq. (^) the master equation for 
quantum Brownian motion derived by Caldeira and Leggett using the oscillator bath model (IJ) (see also JlJ] ) . 
In the case of a Levy stable distribution, the characteristic functional is given by jfjl 



$ L [k(t)] =exp{i J j(t)k(t)dt-J C(t)\k(t)\ a l + i/?A t an^ dt} 



(7) 



where a (0 < a < 2) is the characteristic exponent (or stability index) of the distribution and (3 (—1 < /? < 1) is 
the asymmetry parameter |T(| . We assume as before that ^(t) is proportional to the velocity of the system and take 
C{t) = D constant. This leads to the following quantum master equation for a Levy stable process 



. dp(x,x',t) 
1 Wt 



H(x) - H(x) +2m( x ~ x )(Px -Px')~ 



\x — x\ + i/3tan ~2~i x ~ x ')j p( x : x '>t) 

(8) 



For a = 2 the master equation (|8|) reduces to the Caldeira-Leggett equation (^|). In order to obtain the corresponding 
classical transport equation, we introduce the Wigner transform of the density matrix 



f{q,p,t) 



1 

2~Trh 



dr exp 



ipr 



r r 
p[q+ 2'9 _ 2'* 



(9) 



Applying the Wigner transform to Eq. (g) and keeping only terms in leading order in h, we obtain the equation 

d a f 



df P df , df 7 d 

m=-MT q +u{q) dp: + Mdp- {pf)+lkT 



an d d a - x f 
p tan — 



_d\p\ a ' 2 dpd^" 1 
where we have introduced the Riesz fractional derivative which is defined through its Fourier transform as jL7|Jl 



(10) 



irj^ 2^,' ' / '" M) 



ipr 



(11) 



The fractional equation ( |10[ ) for the distribution function f(q,p,t) describes the complete dynamics of the Brownian 
system in phase-space, for both symmetric and asymmetric Levy stochastic forces. In the latter case, we observe an 
additional contribution to the friction term which may be of relevance for the description of anomalous transport in 
anisotropic media JlS[] . It is worthwhile to note that the fractional character in Eq. ( fl0| ) is carried by the momen- 
tum. For a = 2 we recover the ordinary Klein-Kramers equation. In the limit of high friction, one may exploit the 
rapid relaxation of the momentum distribution to a stationary distribution, to write down a simplified equation for 
the reduced distribution function in configuration space, f(q,t) = J dpf(q,p,t). Employing the systematic expan- 
sion method developed in Ref. |2(i| ], we obtain the following fractional extension of the Smoluchowski equation (the 
derivation of Eq. (|lj) will be sketched below) 



! = ijV( ?) /, + ^ 

at 7 oq 7 



d a f 



d\q[ 



(3 tan • 



Q gat-lf 



2 dqdlq^- 1 



(12) 



We see that the fractional character has been transferred to the position. This observation settles an apparent point of 
confusion in the literature. Indeed, in most of the phenomenological fractional generalizations of transport equations 
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to phase-space, it is unclear whether the fractional derivative should be taken with respect to position or momentum 
or even both (see e.g. the discussion in Refs. [^l|,|2)). Our findings show that for a Levy flight the fractional character 
acquired by the position in the Smoluchowski equation follows from the fractional character of the momentum in the 
Klein-Kramers equation. 

Let us now discuss the connections of the fractional Smoluchowski equation (|l^) to the equations considered in the 
literature. We first begin with the case of a symmetric Levy force. Equation (|l^) with (3 = (which in this form is also 
sometimes called fractional Fokker-Planck equation) has been studied in detail in Ref. Q for the cases of a free flight, 
a particle subjected to a constant force and to a linear Hookean force. Moreover, it is interesting to note that Eq. ( |l2] ) 
has been recently derived in Ref. ]l0| from a generalized master equation for a non-homogeneous random walk. The 
corresponding fractional diffusion equation obtained by setting [/(g) = has also been considered in Ref. . On the 
other hand, for an asymmetric random force, [3 ^ 0, the Smoluchowski equation (|l^) is a generalization to a velocity 
dependent damping force of a fractional diffusion equation recently obtained in Ref. ||] starting from a Langevin-like 
equation. 

In a recent work, Kusnezov et al. [|TJ proposed a fractional Klein-Kramers equation which was obtained as the 
classical limit of what they call a quantum Levy process. Their derivation is based on a microscopic random-matrix 
model for a system coupled to a chaotic environment. These authors showed that for the case of an environment 
with constant average level density (or equivalently with infinite temperature (3t = (fcT) -1 = 0), the reduced density 
matrix displays the behaviour of a free Levy flight. However, the dynamics described by their fractional transport 
equation for finite temperature (3t ^ is unknown. Let us now examine that point. The characteristic functional 
corresponding to their quantum master equation can be written as 



$ QL [k(t)} = exp{i / 7 (t)| BgnfeWlfcCt)! " 1 *- / C(t)\k(t)\ a dt} . (13) 



By comparing expression (|13J) to the characteristic functional of a symmetric Levy flight (|7[), we observe (i) that the 
second terms on the rhs, those describing the fluctuation properties of the stochastic force, are equal, but (ii) that 
the first terms, which are related to the mean, are different. Since the latter are responsible for the dissipation, this 
implies that two expressions are identical only for vanishing friction. Note that this is in agreement with the results 
of Ref. since the limit of vanishing friction precisely corresponds to (3t = 0, as can be easily seen from their 
equation (28). The fractional Kramers equation (27) (with /3t ^ 0) given in Ref. |2l| thus describes a Brownian 
system subjected to a symmetric Levy stable random force, but with a mean friction force that is different from that 
of a Levy flight. It is straightforward to determine the mean value (-FqlW) for the process defined by Eq. (|l3|). It is 
given at a particulary time t — to by 



Sk(t ) 



(14) 



fc=0 



Since a — 1 < 1, the mean force (-Fql(£o)) is divergent (it is finite and equal to 7(^0) only for a — 2). Hence both the 
first and the second moment of the process investigated by Kusnezov et al. are divergent. This has to be contrasted 
with the normal Levy flight, where the mean is finite and equal to (Fi,(to)} = 7(^0) for all values of a. Let us emphasize 
that the random-matrix results are therefore at variance with the continuous time random walk calculation fiofl . The 
two approaches only agree for the special case of a Gaussian distribution a — 2. The divergence of the damping force 
in the random-matrix model can be seen to be related to the Markovian property of the transport equation. 

We now return to the derivation of Eq. ( [l2| ) . For simplicity we will only consider the case of a symmetric probability 
distribution. For large 7, the dynamics of the Klein-Kramers equation ( fTo| ) is dominated by the term which contains 
the operator 

C=-^-p + MkT^ r . (15) 
dp d\p\ a v ' 

We shall look for an approximation to leading order in (7/M) -1 of Eq. ( |Io| ) by using an eigenvalue method |20|| . We 
denote by ip n (p) the eigenfunctions of the operator C and by — n (n = 0, 1, 2, ...) the corresponding eigenvalues. We 
define the following raising and lowering operators 

a+ = -kT^- and a. = MD^ + ^p, (16) 

where the operator obeys — . The two operators a + and a_ satisfy C — — a+a_ and [a_, a+] = 1 and 

we have further the ladder relations a + Lp n {p) = (n + 1) (p n +i(p) and a-(p n (p) = (1 — S n fl) ip n ~i(p). Next we look for 
a solution of Eq. ( [Toj ) in the form 
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f(q,p,t) = f(q,t) <Po(p) 



1_ 

M 



f (1) (q,P,t) + 



M 



f {2 \q,P,t) + .. 



(17) 



and 



df(q,t) 

dt 



= 0<°> + 



M 



9 (D 



M 



a (a) 



/(<?,*) , 



(18) 



where the c?W are linear differential operators which are determined as follows: we substitute Eqs. ( p7| ) and ( p"8| ) into 
Eq. ( |l0| ) and separate the different orders in ('y/M)~ 1 . The integrability condition then yields for the two lowest 
orders 



d<n = o 



and 



a™ = 



kT d 
~M~dq 



D' 



kT 



U'(q) 



(19) 



The operator 9^ is precisely the one appearing in the Smoluchowski equation (|l^). 

To conclude, using influence functional methods, we gave the first derivation of a quantum master equation for 
symmetric and asymmetric Levy flights with viscous damping. By taking the classical limit, we then obtained an 
extension of the Klein-Kramers equation containing fractional derivatives with respect to momentum. In the limit of 
strong damping, this equation was shown to reduce to a fractional Smoluchowski equation with fractional derivatives 
with respect to position. Furthermore, for symmetric Levy stable laws, our results are in agreement with those of 
Kusnezov et al. in the limit of vanishingly small friction. For non-zero friction, we found that the process described by 
their fractional Kramers equation possesses a divergent mean damping force. In the opposite limit of strong friction, 
we recovered the fractional Fokker-Planck equations considered in Refs. |(|[l(|. Finally, for the case of asymmetric 
Levy stable distributions, we gave an extension to a velocity dependent damping force of a fractional diffusion equation 
suggested in Ref. M. 
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